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Abstract 

The isotropic elastic moduli closest to a given anisotropic elasticity tensor are de- 
fined using three definitions of elastic distance, the standard Frobenius (Euclidean) 
norm, the Riemannian distance for tensors, and the log-Euclidean norm. The closest 
moduli are unique for the Riemannian and the log-Euclidean norms, independent of 
whether the difference in stiffness or compliance is considered. Explicit expressions for 
the closest bulk and shear moduli are presented for cubic materials, and an algorithm 
is described for finding them for materials with arbitrary anisotropy. The method is 
illustrated by application to a variety of materials, which are ranked according to their 
distance from isotropy. 

1 Introduction 

The objective here is to answer the question: what is the isotropic material closest to a given 
anisotropic material. In order to attempt an answer one needs a distance or length function 
which measures the difference between the elastic moduli of two materials. The Euclidean 
norm provides a natural definition for distance, and using it one can find the elastic tensor 
of a given symmetry nearest to an anisotropic elastic tensor pp, see also [21 HE IH El El- 
The Euclidean distance function is, however, not invariant under inversion, i.e. considering 
compliance instead of stiffness, and as such does not lead to a unique answer to the question 
posed. To see this, let 5Cijki and SSijki be the elements of the fourth order tensors for the 
differences in elastic stiffness and compliance, respectively. Define the length of a fourth 
order tensor with elements Tijki by (TijkiTij^i) 1 ' 2 . Then it is clear that the length using 
SCijki is not simply related to that of SS^m- 

Recently and separately, Moakher [7] and Arsigny et al. [S] (see also [§]) introduced 
two distance functions for elasticity tensors which are unchanged whether one uses stiffness 
or compliance. The two measures of elastic distances, called the Riemannian distance [Zj 
and the log-Euclidean metric [Hj, each provide a means to unambiguously define the distance 
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between any two elasticity tensors. The focus here is on finding the isotropic material closest 
to a given arbitrarily anisotropic material. 

The distance functions are first reviewed in Section El along with the more common 
Frobenius or Euclidean norm. The theory is developed in terms of matrices, with obvious 
application to tensors. Preliminary results for elastic materials are presented in Section 3, 
where closed-form expressions are derived for the isotropic moduli closest to a given material 
of cubic symmetry. The general problem for materials of arbitrary anisotropy is solved in 
Section 4, and applications to sample materials are described in Section 5. 

2 Matrix distance functions 

We begin with V(n), the vector space of positive definite symmetric matrices in M nxn , the 
space ofnxn real matrices. Recall that a matrix P is symmetric if x T Py = y T Px for 
all x, y in IR n , and positive definite if x T Px > for all nonzero x G M n . The spectral 
decomposition is 

n 

P = ^A iVi vf, (1) 

i=l 

where Aj are the eigenvalues and v, G M n the eigenvectors, which satisfy Aj > 0, v?Vj = 5ij. 
Functions of P can be readily found based on the diagonalized form, in particular, the 
logarithm of a matrix is defined as 



n 

LogP = £)lnA i v < vf . (2) 

i=i 



Three distinct metrics for positive definite symmetric matrices are considered: the conven- 
tional Euclidean or Frobenius metric <1f, the log-Euclidean norm d^ |E] and the Riemannian 
distance d# [7j. Thus, for any pair A, B G V{n) 

d F (A,B) = ||A-B||, (3a) 
d L (A,B) = ||Log(A)-Log(B)||, (3b) 
d*(A,B) = ||Log(A -1/2 BA~ 1/2 )|| , (3c) 

where the norm is defined ||M|| = [t^JVP^M)] 1 / 2 for any M G M nxn . The distance function 
d/j is a consequence of the scalar product 

(Mi, M 2 ) P = tr(p- 1 M 1 p- 1 M 2 ), (4) 

for P G V(n) and symmetric Mi, M 2 G M nxn , and is also related to the exponential map 
|1U[ [7j. The metric d^ is associated with the Lie group on V(n) defined by the following 
multiplication that preserves symmetry and positive definiteness jH] : 

Pi P 2 = exp (Log(Pi) + Log(P 2 )) , Pi, P 2 G V(n). (5) 
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The three distance functions possess the properties expected of a distance function d: (i) it 
is symmetric with respect to its arguments, d(A, B) = d(B, A); (ii) non-negative d(A, B) > 
with equality iff A = B; (iii) it is invariant under a change of basis, d(QAQ T , QBQ T ) = 
d(A, B) for all orthogonal Q G M nxn , QQ T = Q r Q = I; and (iv) it satisfies the triangle 
inequality d(A, C) < d(A, B) + d(B, C) for all A, B, C G V{n). 

The Riemannian and log-Euclidean distances have additional properties not shared with 

d F : 

d Lji? (aA,aB) = d L , R (A,B) , a G R+, (6a) 
d LiR (A\B b ) = |6|d L , fl (A,B), &6l, (6b) 

where &l,r signifies either d^ or d R . Thus di and d R are bi-invariant metrics, i.e. distances 
are invariant under multiplication and inversion. This property makes them consistent and 
unambiguous metrics for elasticity tensors. Moakher [7j introduced another bi-invariant 
distance function, the Kullback-Leibler metric, but it does not satisfy the triangle inequality, 
and we do not consider it here. 

The distance function d^ can be expressed in alternative forms by using the property 
B(LogA)B- 1 = Log (BAB" 1 ), e.g., 

d R (A,B) = [trLog^A^B)] 172 = [tr Log^B^A)] 1/2 , (7) 

or in terms of eigenvalues, using eqs. (0) and (J3cJ), 



dfl(A,B) 



^ (In A 



o 2 



1=1 



1/2 



(8) 



where \, i = 1,2, ... ,n are the eigenvalues of P = A 1/l2 BA 1 / 2 , or equivalently, of the 
matrices A -1 B, B _1 A, AB" 1 , etc. Note that d R also satisfies 

d R (SAS T , SBS T ) = d R {A, B) , V invertible S G M nxn . (9) 



3 Preliminary examples 

The remainder of the paper is concerned with applications to elasticity, with n = 6. 
3.1 Definition of elastic moduli 

6x6 symmetric matrices are used to describe elastic moduli, whether of stiffness or com- 
pliance. The matrix representation is based on Kelvin's observation in 1856 that the 
twenty one coefficients of the elasticity define a quadratic form (the energy) in the six strains, 
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and therefore possess six "principal strains" [TTJ. Although Kelvin did not write the elastic- 
ity tensor explicitly as a symmetric positive definite matrix, the idea has proven useful and 
has been developed extensively, notably by Rychlewski and Mehrabadi and Cowin |13j . 
The notation of Mehrabadi and Cowin ^3] is employed here. Thus, the matrix C G P(6) 
represents the elastic stiffness, and its inverse is the elastic compliance, S, satisfying 



SC = CS = I, where I = diag(l, 1, 1, 1, 1, 1). 
The elements of the elastic stiffness matrix are 

/Cll C12 Ci3 Ci4 Ci 5 C 16 \ ( 

Cl2 C 2 2 C 2 3 C 2 4 C 2 5 C 2 6 c 12 C 22 C 2 3 

Cl3 c 23 c 33 c 34 c 35 c 36 c 13 c 23 c 33 

C14 C24 C34 C44 C45 C 46 
C15 C25 C35 C45 C55 C56 



(10) 



c 



\Cl6 C26 C36 C46 C56 C66/ 



\22Ci6 22C 26 22C 36 2c 



22C M 


25Ci 5 


2^c 16 \ 


2^C 2 4 


25C 25 


2^c 26 


2^c 34 




2^c 36 


2C44 


2c 45 


2c 46 


2c 45 


2c 55 


2c 56 


2C46 


2C56 


2c 66 ^ 



(11) 



where c^-, i, j = 1,2, . . .6 are the coefficients in the Voigt notation. 

Before considering materials of arbitrary anisotropy it is instructive to examine the dis- 
tance functions for materials possessing the simplest type of anisotropy: cubic symmetry. 
Materials of cubic symmetry are described by three independent moduli: hatcn = hatci2 = 
hatc 33 , hatcyi = hatc 23 = hatci 3 , hatc^ = hatc 55 = hatc 6e , with the rest equal to zero. The 
three moduli commonly used are the bulk modulus n and the two distinct shear moduli /i 
and rj, which are related to the matrix elements by 



3k = en + 2c 



12- 



2fi = c 44 , 



2r/ = en - C12. 



(12) 



Isotropic materials have only two independent moduli k, fi, and are of the same form as for 
cubic materials with the restriction en — Ci 2 — c 44 = 0, or equivalently, rj — fi. 

A concise notation is used for isotropic and cubic matrices, based upon Walpole's [13] 
general scheme for performing algebra with elasticity tensors. Define the matrices J, K, L 
and M by 



K : 

M 



I- J, 

K — L, 



J = uu T , where u = (^=, ^, 0, 0, 0) T , (13a) 
L = diag(0, 0, 0, 1, 1, 1) . (13b) 



Note that I and J correspond, respectively, to the fourth order isotropic symmetric tensors 
with components l^i = (SikSji + 5u5jk)/2 and J^i = (l/3)5^5fc/. Elastic moduli of isotropic 
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Ciso(3/t, 2[i) 
C cub (3/t, 2/i, 2rf) 



3nJ + 2fiK, 

3kJ + 2/xL + 2//M, 



k, n > , 
k, /i, r\ > 



(14a) 
(14b) 



The isotropic matrices {J, K} are idempotent and their matrix product vanishes, J 2 = J, 
K 2 = K, JK = KJ = 0. Similarly, it may be checked that the three basis matrices for cubic 
materials {J, L, M} are idempotent and have zero mutual products. The algebra of matrix 
multiplication for isotropic and cubic materials follows from these basic multiplication tables: 





J 


K 


J 


J 





K 





K 





J 


L 


M 


J 


J 








L 





L 


' 


M 








M 



Thus, inverses are S cub = = C cub (^, ^, ^), S iso 



n 3 ^> and products are 



C is ^(3«i, 2/ii)C iso (3fi:2, 2/i 2 ) 
C^uU 3 ^ 2/ii, 2?7 1 )C cub (3K 2 , 2/x 2 , 2?7 2 ) 



K-2 



/'l 
/'2 



(15a) 
(15b) 



K X /ii rji 

Results for isotropic materials follow from those for cubic with r] = \i. For the sake of simplic- 
ity and brevity we therefore focus on properties for cubic materials in the next subsection. 

3.2 Elastic distance for cubic and isotropic materials 

Consider two cubic materials with moduli Ci = C cub (3fti, 2 fix, 2rjx) and C 2 = C cub (3« 2 , 
2/j, 2 , 2?7 2 ). The Euclidean distance function of eq. ()3a|) follows from the above properties 
and the relations trJ = 1, trL = 3, trM = 2. Similarly, the Riemannian and log- Euclidean 
distances follows from the identities 



Log(C 2 ) - Log(Ci) = Log (VC 2 = In (^) J + In (^) L + In (^) M 

Ki \lx Vl 



•1*2- 



>V2- 



Thus, the distances functions are 



d F (Cx, C 2 ) 
d LtR (Ci, C 2 ) 



[(3«i - 3k 2 ) 2 + 3(2//! - 2/i 2 ) 2 + 2(2771 - 2r l2 f] 

K 2x 2 . on ^2, 2 , on ^2x21 1/2 



21 1/2 



[(ln-) 2 + 3(ln 

Kx Hi 



2(ln^) 2 ] 



(16) 



(17a) 
(17b) 



It is clear that d^ and d# are invariant under inversion, di^(Si, S 2 ) = dz^Cx, C 2 ). Note 
that the first identity in ()16j) is a consequence of the fact that Ci and C 2 commute, which 
is not true in general for material symmetries lower than cubic. 
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What is the isotropic material closest to a given cubic material? The answer may be found 
by considering the distance functions between an arbitrary cubic stiffness C cu b(3K, 2/x, 2rj) 
and the isotropic stiffness Cis (3«;*, 2//*). The same question will also be considered for the 
compliances. Minimizing with respect to the isotropic moduli fi* yields 



min (1l,r I ) = min d L>R [ S cub , S 

is' 



min d F [ C cub , C iso (3«;*, 

K*, 



min d F ( C cub , 0^(3^,2/^ 



5 |2M- 277| 



In 



1 

2^ 



1 

2t? 



(18a) 
(18b) 
(18c) 



Denote the values of the closest isotropic moduli by u. l ), (kr, (j,r) for d^, d#, and (ka, Ha) 
or (khj^h) for dp depending as the stiffness (A) or its inverse (H) is used. Thus, 



kl,r,a,h = «, AiL.fl = (A* »7 



3 <n 2 W 5 



2 

5?7 



(19) 



1 3 
/i H 5/i 

Equations (fTSj) and f)19|) show clearly that the "closest" isotropic material using the Frobenius 
metric is ambiguous because it depends on whether one uses stiffness or compliance. Each 
gives a different isotropic material since fiH < Hl,r < Ha for /i - ^ 0. The Riemannian 
and log-Euclidean metrics gives the same unique "closest" isotropic material, regardless of 
whether the stiffness or the compliance are used. The fact that they agree is particular to 
the case of cubic symmetry, as noted above, and is not true in general. 

In summary, the closest isotropic material to a given cubic material, in the sense of d# 
and d^, is defined by moduli kr = kl = |(cn + 2ci 2 ) and hr — Hl = \ [c| 4 (en — Ci 2 ) 2 ] 1 
and the distance from isotropy is di : R = 
to materials of arbitrary anisotropy next. 



6/5 



C44 



. These results will be generalized 



4 Closest isotropic moduli 

We now turn to the more general question of finding the isotropic material closest to a given 
anisotropic material characterized by C or its inverse S. The solution using the Euclidean 
metric is relatively simple, and is considered first. 

4.1 Minimum Frobenius distances 

The closest isotropic elastic moduli are assumed to be of general isotropic form Ci SO (3K, 2/x), 
see eq. ()14j) . The bulk and shear moduli are found by minimizing di?(Ci SO , C), which implies 



3KtrJ = trJC, 2/itrK = trKC 



(20) 
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Using suffix A to indicate that the minimization is in the arithmetic sense (in line with [7j), 

9k a = en + c 22 + c 33 + 2(c 23 + c 3 i + C12) , (21a) 
30fi A = 2(cn + c 22 + c 33 - c 23 - c 3 i - C12) + 3(c 44 + c 55 + c 66 ) , (21b) 

which are well known, e.g. Similarly, the closest isotropic elastic compliance can be 

determined by minimizing dp(C^, C _1 ). Denoting the isotropic moduli with the suffix H 
for harmonic, 

\/k h = sn + s 22 + s 33 + 2(s 23 + s 3 i + S12) , (22a) 
15/ (2/Ih) = 2(s n + s 22 + s 33 - S23 - S31 - S12) + 3(s 44 + S55 + s m ) ■ (22b) 

The Euclidean distance does not provide a unique closest isotropic material, although the 
values in eqs. (|2*Tf and (j22J) are sometimes considered as bounds. Equations (|2T)j) and (|2"Tf 
also agree with the special case discussed above for cubic materials, eq. (fT9|) . 

4.2 Minimum log-Euclidean distance 

The isotropic elasticity Ci SO (3/€L, 2/i^) is found using the same methods as above by replacing 
C iso and C with Log(C iso ) and Log(C), respectively. Thus, 

log(3K L ) = tr JLog(C), 5 log(2/i L ) = tr KLog(C) . (23) 

Adding the two equations and using J + K = I, implies the identity 

det(C iso ) = det(C). (24) 

Thus, we have explicit formulae for the closest moduli, 

k l = § exp (tr JLog(C)) , fi L = | exp (| tr KLog(C)) . (25) 

4.3 The minimum Riemannian distance 

We look for moduli of the form 

C 1so (3k r , 2/i r ) = 3k r J + 2jj, R K , (26) 

which minimize dij(C; so , C). Differentiating the expression 

d^(C iso , C) = tr [LogCCr^C) Log(Cr o 1 C)] , (27) 

with respect to n R and /ir separately, and using Log(A)' = A" 1 A' for A G V(n), implies 
respectively 

tr [C^JLog(C^C)] =0, tr [c£KLog(C£C)] =0. (28) 
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Further simplification yields 

tr pLog(CiC)] =0, tr [K Log(Cj^C)] =0. (29) 

These conditions, which are necessary for a minimum, can be simplified as follows. Define 
the eigenvalues and associated eigenvectors by the diagonalization 

n 

C^CC^^Wf. (30) 
t=i 

Adding the two conditions (J2TJJ) using the identity I = J + K, along with the expression (0) 
for the logarithm of a matrix, yields 

n 

n^ =i - (3i) 
i=i 

A second condition follows by direct substitution from (130(1 into (|29|)i . giving 

n 

H Af* = 1, c^vfjv,, z = l,2,. ..n. (32) 
t=i 

Note that < a,, < 1 and ccj form a partition of unity, 

n 

<*i = 1 ■ (33) 

i=l 

This follows from the representation J = uu T where the unit 6— vector u is defined in eq. 
(j!3a|) . Thus, the minimal isotropic moduli are found by satisfying the two simultaneous 
equations (|31jl and fl32(l . We now show how the first of these two conditions can be met, 
leaving one condition to satisfy. 
Let 

C iso = 3k r (J + p~ 2 K) , (34) 

where p > is defined by 



3k r 1 + v R 



and vr is the Poisson's ratio of the minimizer. The reason this form for Ci SO is chosen is so 
that C^y 2 = (3K R y^ 2 (j + pK). Hence, the eigenvalues of §^ are of the form 

K = -5 , (36) 

3k r 

where the normalized eigenvectors \ = Xi(p), and the (unchanged) eigenvectors Vj, i — 
1, 2, . . . , n — 6 are defined by 

3 Ki? Cr o 1/2 CCr o 1/2 = ( J + P K) C{3 + P K)=J2 \ Vivf . (37) 

8=1 
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Turning to the first condition, (|31|). it is automatically satisfied if the bulk modulus is 
given by 

/ n _ \ i/« 

3k r =(Y[\A . (38) 
^ i=i ' 

It remains to determine p from the second stationary condition, eq. ()32j). which can be 
expressed in terms of the modified eigenvalues as 

n 

J! A^- 1/n) = 1 . (39) 

i=l 

Equation ((Hi"") involves the eigenvectors v through the inner products c^. However, on van- 
ishes identically for eigenvectors of deviatoric form - in fact the definition of a deviatoric 
eigenvector is = [12]. Conversely, «j = 1 for purely dilatational eigenvectors [12], i.e., 
eigenvectors parallel to u of eq. (J13a() . 

The solution to eq. (J39)) may be found numerically by searching for the zero over the 
permissible range for the Poisson's ratio: — 1 < v R < 1/2. The minimizing moduli k r and 
p R then follow from eqs. (J3~8*|) and ((3l"") . or more directly, 

3k r = p 5 / 3 (det C) 1/6 , 2p R = p- 1/3 (det C) 1/6 , (40) 

and the minimal distance between Ci SO and C is given by 

1 



di?,(Ci SO , C) — — 
n 



$:in 2 (Ar-n h\ 



1/2 

(n = 6). (41) 

L i=i j=i 

We next demonstrate the application of the above procedure to the case of a given 
elasticity matrix of cubic symmetry. 



4.4 Example: cubic materials 

Substituting the assumed form C = C cu b from eq. (j!4b|) into the explicit formulae of eq. 
(125(1 for the closest moduli in the log-Euclidean sense, it is a straightforward matter to show 
that these reproduce the results determined directly, in eq. (|19(1. Regarding the closest 
moduli using the Riemannian distance, the matrix in eq. (|37(1 follows by using the algebra 
for cubic matrices, 

(J + pK) C (p J + pK) = 3k J + 2/ip 2 L + 2r]p 2 M . (42) 

Thus, Ai = 3k, \ 2 = A3 = A4 = 2pp 2 , X5 = \e = 2np 2 , and the eigenvectors are either pure 
dilatational (a% = 1) or deviatoric (aij = 0, i = 2, 3, . . . , 6). Therefore, eq. ([3~9~|l becomes 

(3k) 5/6 {2pY l/2 (2r ? )" 1/3 p- 5 / 3 = 1. (43) 

Solving for the intermediate variable p, and evaluating p R and k r from eqs. ([3~5|l and ((3*3(1 
respectively, gives k r = k and p/j = (p 3 // 2 ) 1 / 5 , again in agreement with eq. ("B"() . 
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Table 1 lists the computed distance from isotropy of various anisotropic materials, using 
data from Musgrave [TSJ. Materials of cubic (cub), hexagonal (hex), tetragonal (hex) and 
orthotropic (ort) symmetry are considered. In each case the moduli of the closest isotropic 
material were found using the algorithm described above. The resulting bulk modulus k,r 
and Poisson's ratio Vr are tabulated. 

Table 1 ranks the materials in terms of the Riemannian distance d^ of the original 
anisotropic moduli from the closest isotropic material. The second column of numbers lists 
the distance between the closest isotropic materials found using the Riemannian and log- 
Euclidean norms. That is, 

d LR = d L , R (C iso (3/^,2/i R ), 6^(3^,2^)) = [(In ^) 2 + 5(ln ^) 2 ] 1/2 , (44) 

\Xr 

which is identically zero for cubic materials. The arithmetic (ka, Ha) an d harmonic (kh, Hh) 
moduli which minimize the Euclidean distances were also computed, and the Riemannian 
distance between these two is denoted da a- The distances d^ and dRH are the distances 
between the closest isotropic material (kr, hr) and the arithmetic and harmonic isotropic 
approximants, respectively. All distances listed in Table 1 are based on the Riemannian 
metric. 

Note that the distance between the closest materials using d.R and d^ is less than 0.05 
except for the extremely anisotropic spruce. In order to gain some appreciation for the 
magnitude of the nondimensional distances in Table 1, consider the distance of any P £ V(n) 
from a multiple of itself: 

d R (P, aP) = d L (P, aP) = ^ |loga| , a £ R+ . (45) 

Small values of the elastic distance can be identified with values of a close to unity, specifically 

a = 1 ± ^d^R + 0{d 2 L)R ) « 1 ± 0.4 d L , R . (46) 

Note that the distance duA between the arithmetic and harmonic approximations is 
generally less than the distance from isotropy d#. This is more so for those materials that 
are closer to isotropy - at the top of Table 1. As the material gets further from isotropy - 
the lower half of Table 1 - the magnitude of duA relative to dR grows as the latter increases. 
The two distances are of comparable magnitude for the highly anisotropic materials at the 
very bottom of the table, such as oak and spruce. 

As a numerical check on the computations, the triangle inequality 



dnA < dRA + dRH, 



(47) 
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was confirmed for each material in Table 1. Since the three vertices of the triangle are 
isotropic materials, the inequality may be written, using (|17b[) . as 

[(In ^-f + 5(ln ^f] 1/2 < [(In + 5(ln ^) 2 ] 1/2 + [(In + 5(ln ^f] 1/2 . (48) 

For cubic materials k,a = k h = K Ri and consequently the triangle is a straight line: 

&HA = d^ + dun for cubic materials only. (49) 

The quantity (cIra + ^rh — &ha)I&ha wa s found to be very small for all the cases 
considered (and numerically zero for the cubic examples), less than 10~ 3 for all materials 
considered except barium titanate (1.2x 10~ 3 ) and spruce (2.8x 10~ 3 ). The "triangle" is 
almost flat, indicating that the closest moduli (kr, /ir) are in some sense optimally centered 
between the arithmetic and harmonic approximations. Note however, that kr, fiR are not 
equal to the Riemannian mean [Z| of the arithmetic and harmonic approximations, denoted 
as Hah-, Hah- The Riemannian mean of two elasticity matrices Ci and C2 is Ci(C]~ 1 C2) 1 ^ 2 
[7j, and consequently the means of the arithmetic and harmonic moduli are k^h — ( k a k hY^ 2 , 
Hah = (ha^h) 1 ^ 2 - By considering the case of cubic materials, for which all these quantities 
have explicit expressions, it may be shown that (fiR — hah)(v — ji) > for r\ — /1 ^ 0. 



Conclusions 

We have presented a method for finding the isotropic elastic moduli closest to a given material 
of arbitrary symmetry based on three different metrics. Unlike the Frobenius (Euclidean) 
norm, the Riemmanian and log-Euclidean metrics provide unique isotropic moduli. The 
values obtained according to these two metrics are identical if the comparison medium has 
cubic symmetry, and are otherwise relatively close. The procedures developed here for 
finding the closest isotropic moduli can be generalized to find the closest material of lower 
symmetry. The solution for cubic symmetry with the cube axes given is presented in the 
Appendix, and other, lower, symmetries will be considered elsewhere. Another generalization 
of the present problem is that of determining the closest material of cubic or lower symmetry 
where the symmetry axes are unrestrained. These and other challenging questions make this 
an interesting topic for some time to come. 
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Appendix 

A. The closest cubic material 

The cubic stiffness (compliance) closest to C (S) in the Euclidean metric dp has moduli k^, 
Pa, Va (kh, PiT and t]h), where ka and Kh are given by ()21a|) and (|22a|) . and 

6/iA = c 44 + c 55 + c 66 , 6p A = cu + c 22 + c 33 - c 23 - c 3 i - C12, (A.la) 
3/ (2p#) = -§44 + S55 + s 6 6, 3/ (277h) = sn + s 22 + s 33 - S23 - s 3 i - S12 . (A. lb) 

For the log- Euclidean norm, we find, using the method for deriving eq. (|25|). 

k l = I exp (tr JLog(C)), p L = § exp (-trLLog(C)), pl = | exp (| tr MLog(C)) . 

(A.2) 

Note the identity, similar to (J23|) . 

det(C cub ) = det(C). (A.3) 

For the Riemannian norm d# we find that the closest cubic material C cu b of the form 
()14b|) is determined by three equations: 

n n n 

n A * =i ' n A " =i ' rK =i ( a - 4 ) 

i=l i=l i=l 

where 

an = vf Jvi, A = vf Lvi, i = 1, 2, . . . n, (A. 5) 

■ — ■ ]/2 """^ 1/2 

and {Aj, Vj} are the eigenvalues and eigenvectors of C cub CC cub . The parameters a, satisfy 
the same properties as before, including the fact that they sum to unity. Since {vj} form 
an orthonormal basis, it follows that Y^=i A = dimL = 3. Furthermore, fa — if the 
eigenvector is dilatational. The three equations (jA.4|) may be reduced to two by assuming 
the unknown moduli are of the form C cu b = 3kr (j + p]~ 2 L + 2 M) . Define the modified 
eigenvalues \ = Xi(pi, P2) to be the eigenvalues of (J + piL + P2M) C ( J + piL + P2M) , then 
k r is given by the formula while pi, pi solve the simultaneous equations 

n n 

n A 4 K - i/n) =1, n Af- i/n) = 1 . (a.6) 

i=l i=l 
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Table 1: Distance from isotropy for some materials - data from ^H]- k r units 10 10 N/m 2 



Material 


bymm 


d R 


100d L/? 


d 


RA 


d 


RH 


d 


HA 








magnesium 


1 

hex 





,18 


0.00 


0. 


01 


0, 


01 


0, 


.02 





,29 


3 


,53 


diamond 


cub 





.21 





0, 


,01 


0, 


01 


0. 


02 





,07 


44. 


,20 


aluminum 


cub 





.21 





0, 


,01 


0, 


01 


0. 


02 





,35 


7 


.69 


1 11* 

beryllium 


hex 





.22 


0.01 


0, 


,01 


0, 


01 


0. 


02 





,05 


11 


.44 


sodium fluoride 


cub 





.29 





0, 


02 


0, 


02 


0. 


,04 





,24 


4 


.86 


ice (H 2 0) 257°K 


hex 





.31 


0.00 


0, 


02 


0, 


02 


0. 


04 





,33 





.89 


/3-quartz (Si0 2 ) 


hex 





.35 


0.02 


0, 


03 


0, 


03 


0. 


05 





,21 


5 


.64 


beryllium 


hex 





,37 


0.23 


0, 


03 


0, 


03 


0. 


06 





,26 


14 


,41 


caesium iodide 


cub 





,37 





0, 


03 


0, 


03 


0, 


06 





,27 


1 


,29 


sodium chloride 


cub 





,40 





0, 


04 


0, 


03 


0. 


07 





,25 


2 


,45 


sodium iodide 


cub 





,43 





0, 


04 


0, 


04 


0, 


.08 





,25 


1 


,46 


i • 1 * 1 

sodium bromide 


cub 





,44 





0, 


04 


0, 


.04 


0. 


09 





,25 


1 


,94 


1 * i 

caesium bromide 


cub 





,45 





0, 


05 


0, 


.04 


0. 


09 





,27 


1 


,59 


• 1 • 

silicon 


cub 





,49 





0, 


05 


0, 


.05 


0. 


11 





,22 


9 


,78 


cobalt 


hex 





,51 


0.00 


0, 


07 


0, 


,05 


0. 


12 





,31 


19 


,03 


• 1 1 ■ 1 

silver bromide 


cub 





.52 





0, 


06 


0, 


06 


0. 


12 





,40 


i 


.06 


germanium 


cub 





.56 





0, 


,07 


0, 


07 


0. 


14 





,21 


7 


.52 


caesium chloride 


cub 





.58 





0, 


,08 


0, 


07 


0. 


15 





,27 


1 


.83 


gallium antimonide 


cub 





,64 





0, 


09 


0, 


10 


0. 


18 





,25 


5 


,64 


a-uranium 


ort 





,68 


0.37 


0, 


10 


0, 


10 


0, 


20 





,20 


11 


,28 


silver chloride 


cub 





,70 





0, 


11 


0, 


10 


0, 


.22 





,41 


4 


,42 


apatite 


hex 





,72 


0.11 


0, 


10 


0, 


13 


0, 


.22 





,21 


8 


,43 


1 * j • • i 

indium antimonide 


cub 





,75 





0, 


12 


0, 


13 


0, 


25 





,29 


4 


,69 


potassium fluoride 


cub 





,75 





0, 


13 


0, 


12 


0, 


25 





,28 


3 


,19 


benzophenone 


ort 





,85 


1.92 


0, 


15 


0, 


14 


0, 


29 





,30 


5 


,14 


zircon 


tet 





,98 


0.74 


0, 


21 


0, 


18 


0, 


39 





,13 


1 


,99 


sulphur 


ort 





.98 


4.13 


0, 


20 


0, 


18 


0. 


39 





,34 


1 


.88 


iron 


cub 





.99 





0, 


20 


0, 


23 


0. 


43 





,30 


17 


.05 


* 1 1 

nickel 


cub 


1 


,01 





0, 


21 


0, 


23 


0. 


44 





,29 


18 


.04 


cadmium 


hex 


1 


.04 


3.43 


0, 


20 


0, 


.24 


0. 


,44 





,30 


5 


.40 


rutile (Ti0 2 ) 


tet 


1 


,07 


0.79 


0, 


21 


0, 


.28 


0, 


,49 





,27 


21 


.49 


potassium chloride 


cub 


1 


,08 





0, 


27 


0, 


.24 


0, 


50 





,28 


1 


,78 


barium titanate 


tet 


1 


,13 


3.20 


0, 


,26 


0, 


.27 


0. 


52 





,36 


17 


,67 


potassium bromide 


cub 


1 


,14 





0, 


30 


0, 


26 


0, 


56 





,29 


1 


,58 


gold 


cub 


1 


,16 





0, 


.27 


0, 


31 


0. 


.58 





,42 


17 


,28 


T— V 111 1 1 

Kochelie salt 


ort 


1 


,17 


0.97 


0, 


.24 


0, 


34 


0, 


59 





,31 


1 


,97 


zinc 


hex 


1 


"1 o 

.18 


2.58 


0, 


,24 


0, 


.34 


0. 


,57 





,24 


6 


.61 


white tin 


tet 


1 


,18 


0.04 


0, 


,24 


0, 


■~> n 

.38 


0. 


,62 





,35 


5 


.50 


ammon. dihyd. phos. 


tet 


1 


,19 


0.95 


0, 


36 


0, 


25 


0. 


61 





,33 


2 


,70 


silver 


cub 


1 


,21 





0, 


29 


0, 


33 


0. 


63 





,37 


10 


.36 


potassium iodide 


cub 


1 


.25 





0, 


36 


0, 


31 


0. 


67 





,30 


1 


.20 


copper 


cub 


1 


,28 





0, 


32 


0, 


37 


0. 


70 





,35 


13 


,71 


potass, dihyd. phos. 


tet 


1 


,34 


0.01 


0, 


.40 


0, 


38 


0, 


78 





,26 


2 


,67 


a-brass 


cub 


1 


,46 





0, 


,41 


0, 


.48 


0, 


90 





,34 


11 


.96 


indium 


tet 


1 


,57 


0.01 


0, 


.50 


0, 


,54 


1, 


04 





,44 


4 


,16 


oak 


ort 


2 


,30 


1.75 


0, 


96 


1 


09 


2. 


,05 





,08 





,17 


/3-brass 


cub 


2 


,34 





0, 


.94 


1, 


,19 


2. 


13 





.36 


11 


,62 


spruce 


ort 


5 


.66 


59.5 


7, 


16 


3 


33 


1C 


i.45 





.23 





.09 



